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Energy  Absorption  of  Graphite/Epoxy  Plates 
Using  Hopkinson  Bar  Impact 

PIYUSH  K.  DUTi  A.  DAVID  HUI  AND  MAGNA  R.  ALTAMIRANO 


INTRODUCTION 

Most  literature  dealing  with  energy  absorption  in 
laminated  composites  describes  the  use  of  a  dropweight 
or  pendulum  impact  test  apparatus.  In  these  techniques, 
energy  absorption  is  determined  by  subtracting  the  re¬ 
sidual  kinetic  energy  from  the  initial  energy.  Thus  the 
effects  of  stress  wave  propagation,  which  is  a  source  of 
damage  initiation,  cannot  be  examined  using  these  test¬ 
ing  techniques.  Further,  the  laminate’s  energy  absorp¬ 
tion  during  the  penetration  process  and  the  projectile’s 
velocity,  contact  force,  and  duration  of  impact  are  diffi¬ 
cult  to  measure.  Hopkinson  bar  testing  eliminates  these 
restrictions  and  allows  a  thorough  examination  of  the 
projectile  impact  process. 

The  effort  described  in  this  report  involved  an  ana¬ 
lytical  and  experimental  study  of  energy  absorption 
within  graphite/epoxy  plates  impacted  by  a  hemispheri¬ 
cal-nosed  impactor  attached  to  a  split  Hopkinson  bar  test 
apparatus.  The  work  is  novel  in  that  the  plate’s  energy 
absorption  was  accurately  measured  using  the  recorded 
incident,  reflected,  and  transmitted  stress  waves.  The 
Hopkinson  bar  test  apparatus  allowed  accurate  determi¬ 
nation  of  velocity,  force,  and  energy  absorption  informa¬ 
tion  during  the  entire  impact  duration.  Thus  the  setup 
allowed  “controlled”  testing,  capable  of  determining  the 
energy  absorption  of  impacted  composite  specimens. 
The  only  1  imitation  was  that  the  impactor  was  forced  into 
the  laminate  at  relatively  low  velocities. 


MECHANICS  OF 
WAVE  PROPAGATION 

The  theoretical  wave  mechanics  involved  in  the  im¬ 
pact  process  are  discussed  based  on  elementary  prin¬ 
ciples  of  unidirectional  stress  wave  propagation.  The 


duration  of  impact  produced  by  the  striker  bar  on  the 
Hopkinson  bar  apparatus  persists  for  a  few  hundred 
microseconds,  and  it  takes  a  small  but  finite  time  for  the 
effects  of  the  impact  force  (stress  wave)  to  enter  the 
graphite  epoxy  test  laminate  at  the  opposite  end  of  the 
bar.  Upon  entering  the  laminate,  a  portion  of  the  incident 
stress  wave  is  propagated  into  the  transmitter  bar  and 
another  portion  is  reflected  back  into  the  incident  bar. 
Incident,  reflected,  and  transmitted  stress  waves  are 
measured  by  strain  gauge  instruments  described  else¬ 
where  (Dutta  et  al,  1987). 

In  the  current  study,  the  laminate  specimens  con¬ 
sisted  of  30  plies  of  AS4/3502  graphite  epoxy  with  the 
following  stacking  sequence: 

[(±  45/02)2/90/+ 45 /02/±  45]^ 

The  indentor  was  driven  into  the  laminate  by  the  force 
of  the  incident  bar’s  stress  wave.  Figure  1  shows  a 
schematic  of  the  Hopkinson  bar,  the  impactor,  and  the 
laminate  panel  clamping  assembly.  Figure  2  is  a  photo¬ 
graph  of  the  setup.  The  transmitter  bar  is  free  to  slide 
forward  or  backward.  As  the  striker  hits  the  incident  bar 
and  drives  the  indentor  into  the  laminate  plate,  the  plate 
bends  and  behaves  like  a  spring-mass  system.  Plate 
bending  was  recorded  from  strain  gauges  mounted  on 
the  rear  face  of  the  laminate. 

The  wave  mechanics  within  this  system  was  analyzed 
to  determine  the  force-displacement  relationship  be¬ 
tween  the  impactor  and  the  composite  laminate.  Energy 
absorbed  during  the  process  was  also  determined  as  a 
function  of  time.  The  velocity,  force,  and  energy  ab¬ 
sorbed  were  calculated  from  incident,  reflected,  and 
transmitted  waveforms. 

An  impact  of  bar  1  (the  striker  bar)  into  bar  2  (the 
incident  bar)  produces  a  compressive  force  at  interface 
A.  According  to  elementary  stress  wave  propagation 


theory  (Kolsky  1963,  Ehitta  et  al.  1987),  stress  a  and 
force  P  generated  by  the  impact  are 

a=(p/g)cV'  (1) 

and 

/'  =  (>l)(p/g){cl/)  (2) 

where  p  is  the  density  of  the  incident  bar,  g  is  the 
acceleration  due  to  gravity,  c  is  the  speed  of  sound  in  the 
bar,  V  is  the  particle  velocity  in  bar  2,  and  A  is  the  cross- 
sectional  area.  The  stress  wave  propagates  along  bar  2  to 
the  opposite  end  where  it  encounters  a  discontinuity  of 
its  impedance  (a  product  of  density  and  stress  wave 
velocity,  pc)  at  interface  B.  In  the  experimental  setup, 
the  impactor  and  Irminate  respond  to  this  incoming 
stress  pulse.  As  a  result  of  the  discontinuity,  a  portion  of 
the  stress  pulse  energy  (£|)  is  reflected  (£g)  and  a  portion 
is  transmitted  (£j)  through  the  spring-mass  system. 
Most  of  the  remaining  energy  goes  into  the  formation  of 
plastic  deformation  of  the  plate.  Some  energy  is  also  lost 
in  overcoming  frictional  forces  in  various  joints  and 
through  noise  and  vibration  of  the  entire  system.  Net 
energy  lost  in  the  impact  can  be  calculated  by  subtracting 
the  energy  of  the  transmitted  and  reflected  waves  from 
that  of  the  incident  wave: 

£a  =  £i  -  £r  -  (3) 


The  force  of  penetration  at  any  time  t  can  be  calcu¬ 
lated  more  directly  by  considering  the  particle  velocity 
in  the  following  equation: 


V  =du/ 


(4) 


where  u  is  the  axial  displacement  of  the  bar  and  is  given 
by 


I 


Under  the  incident  stress  wave  the  penetrator 
end  of  bar  2  will  be  subjected  to  a  displacement  given 
by 


t 


t  =  0 


(6) 


Due  to  the  change  of  impedance,  a  portion  of  the  stress 
wave  is  reflected,  which  produces  a  corresponding  dis¬ 
placement  of  the  indentor  in  the  opposite  direction;  this 
displacement  is 


The  procedure  for  computing  energy  from  the  stress  vs 
time  waveform  is  given  in  Appendix  A. 
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Figure  2.  Photograph  of  the  setup  showing 
the  hemispherical  indentor. 
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The  resultant  displacement  is  given  by 

(8) 

1=0 

Force  on  the  penetrator  [/=■(/)]  is  the  algebraic  sum  of 
the  incident  force  AOjCr)  and  the  reflected  force  Aa/r): 

t 

F(r)  =  A  {[ai(r) +  ar(/)])^f/  (9) 

/=0 

From  eq  8  and  9,  the  force-displacement  curve  forthe 
indentation  can  be  plotted.  The  area  under  the  curve  will 
yield  the  energy  absorbed  as  a  function  of  time  through¬ 
out  the  impact. 

Energy  absorbed  by  the  graphite  epoxy  plate  can  be 
computed  from  the  following  energy  balance  equation 
(Greszczuk  1982,  Shivakumaretal.  1984): 

Ea  =  Ec  +  E\j  +  Em  +  Ed  +  Ef  (10) 

where  =  the  energy  expended  in  Hertzian  contact 
indentation, 

£(,  =  the  energy  due  to  plate  bending, 

E^  =  the  energy  due  to  the  stretching  of  the 
middle  plate  surface, 

Ejj  =  the  energy  due  to  permanent  plate  damage, 
and 

Ef  =  the  energy  due  to  friction,  heating,  etc. 

All  of  these  energies  are  recoverable,  except  the 
damage  and  frictional  energy.  Plate  damage  may  include 
any  or  all  of  the  following  modes  of  failure:  matrix 
cracking,  delamination,  fiber  breakage,  and  spalling. 

The  bending  energy  of  a  circular  plate  clamped  along 
its  edge  under  a  concentrated  load  P  is  (Roark  and  Young 
1975,  Greszczuk  1982) 

Eb  =(P^)(3aV2)(l  -v2)(4nE,/i2j-'  (H) 

where  a  is  the  radius  of  the  circular  plate,  h  is  the 
thickness  of  the  plate,  and  E^  and  are  the  Young’s 
modulus  and  Poisson’s  ratio,  respectively,  of  the  plate 


material  in  the  in-plane  radial  direction.  Since  the  deflec¬ 
tion  was  measured  to  be  less  than  20%  of  the  laminate 
thickness,  the  membrane  energy  (E^^,)  was  neglected. 

Although  Hertzian  contact  law  was  established  for 
local  static  indentation  of  two  elastic  bodies,  it  is  com¬ 
monly  applied  to  impact  situations  where  permanent 
plastic  deformations  occur  (Goldsmith  1960).  The  Hert¬ 
zian  contact  energy  of  a  spherical  indentor  with  radius  R 
impacting  a  flat  plate  is  (Greszczuk  1982) 

£:j.  =  j/>5/3j  (4/5)  (371/16)2/3 

[{k,  +  k2)l[RI2y‘^Y'^ 

(12) 

where 

=(l  -vi)/(7i£s) 

(13) 

*2  =  (l  -v?)/(7tEz) 

(14) 

where  E^  and  are  the  Young’s  modulus  and  Poisson’s 
ratioof  the  steel  indentor,  respectively, E^  is  the  laminate’s 
modulus  in  the  thickness  direction,  and  is  the  laminate ’s 
Poisson’s  ratio. 

Assuming  point  contact  between  the  indentor  and 
plate,  the  deflection  can  be  computed  from  the  in-plane 
stress  a  (Roark  and  Young  1975): 

a/j^/6  =  [E/{4Jt)]  (l  +  v) 

{ln[a/(0.325A)]}  (15) 

Further,  deflection  y  at  the  center  of  the  circular  plate  due 
to  a  concentrated  load  P  is 

y  =  (Pa2/7t)(3/4)[(l  -v2)/(e/i3)]  (16) 

Measured  strain  can  be  converted  to  stress  by  multi¬ 
plying  the  strain  by  the  in-plane  Young’s  modulus,  so  a 
is  a  known  quantity.  Eliminating  the  applied  loadP  from 
the  above  two  equations,  the  deflection  can  be  obtained 
in  the  form 

a^o  (l  -  v) 

y  = - ^ -  (17) 

(2hE)  In  [a/(0.325>i)] 

Upon  stress  wave  arrival,  the  indentor  accelerates 
into  the  plate  and  creates  a  compressive  stress  wave  that 
propagates  through  the  thickness.  Stress  wave  velocity 
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Figure  3.  Stress  pulse  reflections  and  trans¬ 
missions  through  the  interface.  Figure  4.  Interface  velocity  and  force. 


in  the  thickness  direction  of  the  plate  is  approximately 
(McMaster  1963) 

Q  =  [fz/(p/^)'^V(v)]  (18) 

/•(V)  ={(l-v)/[(l  + v)(1-2v)]|'^2 

where  is  Young’s  modulus  in  the  thickness  direction 
and  V  is  Poisson’s  ratio(for  n  =  0.31,^v)  =  1.17732). 
For  a  quasi- isotropic  graphite  epoxy  laminate,  where 
=  5 1 . 1 6  GPa  (7.42  Msi),  =  1 1 .86  GPa  ( 1 .72  Msi), 
and  p  =  1 522.4  kg  m“^  (0.055  lb-mass  in.~^),  wave  speed 
c.^  -  3 .043  km  s" '  (0. 1 1 98  X 10^  in.  s“ ' ).  This  wave  will 
propagate  with  speed  c\  as  a  compressive  wave  and  will 
eventually  reflect  off  the  rear  (free)  surface  as  a  tensile 
wave  with  the  same  velocity.  When  the  wave  arrives  at 
the  rear  surface,  the  plate  moves  away  from  the  indentor 


at  velocity  Yj  while  the  indentor  continues  to  move  into 
the  plate  with  a  velocity  Vj  (see  Fig.  3). 

During  the  duration  of  the  stress  wave,  three  condi¬ 
tions  are  possible.  First,  when  V2<V^ ,  the  indentor  accel¬ 
erates  into  the  plate  (positive  contact  force)  with  a  force 
vs  time  relationship  as  shown  in  Figure  4.  Second,  when 
V2  -  P|,  both  the  indentor  and  the  plate  move  in  the 
incident  direction  with  no  change  in  contact  force.  Third, 
when  V2  >V,,  loss  of  contact  occurs  between  the  in¬ 
dentor  and  plate,  and  a  tensile  force  develops  in  the 
incident  bar.  Both  the  indentor  and  the  plate  move  for¬ 
ward  monotonically  even  though  loss  of  contact  occurs. 

DISCUSSION  OF  RESULTS 

Figure  5  is  a  graph  of  energy  absorption  vs  displace¬ 
ment.  The  energy  absorption  of  the  plate  (curve  A)  is 
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Figure  5.  Impact  energy  absorp¬ 
tion  results. 
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Time  (ns) 

Figure  6.  Results  of  force,  velocity,  and  strain  history  measurements. 


computed  from  the  force  vs  displacement  of  the  indentor. 
Energy  absorbed  in  the  system,  consisting  of  the  plate 
and  clamping  fixture  shown  in  curve  B,  is  computed 
from  the  incident,  reflected,  and  transmitted  waves. 
Appropriate  formulae  are  given  in  Appendix  A.  The 
small  difference  between  these  two  energy  absorptions 
is  the  energy  loss  in  the  fixture  (£  ). 

Indentor  force  vs  time  (curve  C  of  Figure  6)  follows 
the  model  of  Figure  4.  There  are  five  stages  involved  in 
the  penetration  process: 

•  In  the  initial  loading  stage  (0  <  /  <  r,),  the 
indentor  accelerates  into  the  plate,  as  can  be 
seen  from  the  increasing  slope  of  the  velocity 
trace  (curve  D). 

•  The  second  stage  (f,<  r  <  t^)  is  the  initial 
unloading  stage  where  the  indentor  force  de¬ 
creases  rapidly  and  the  corresponding  veloc¬ 
ity  increases  at  a  slower  rate.  The  drop  in  the 
i  ndentor  force  is  attributed  to  damage  either  in 
the  contact  region  or  near  the  rear  surface  due 
to  a  reflected  tensile  wave.  The  latter  may 
cause  the  plate  to  move  away  from  the  indentor, 
thus  reducing  the  contact  force  (which  will  be 
more  evident  in  the  next  stage). 

•  The  third  stage  is  the  constant  velocity  stage 
(/2  <  f  <  ^3)  in  which  there  is  very  little  contact 
force.  Referring  to  eq  9,  note  that  at  this  stage 
aj(f)  =  Oj^(r),  indicating  that  there  is  no  inter¬ 
facial  contact  force. 

•  In  the  fourth  stage  (tj  <  r  <  r^),  the  effects  of 
stage  3  become  more  pronounced,  and  sepa¬ 
ration  between  the  indentor  and  plate  is  pos¬ 
sible.  Similar  loss  of  contact  behavior  has 
been  reported  by  Wu  and  Springer  ( 1 988)  and 
Sun  and  Chen  (1985).  At  this  stage,  the  ap¬ 
plied  incident  stress  amplitude  begins  to  di¬ 


minish  but,  because  of  the  gap  between  the 
plate  and  the  indentor,  the  indentor  experi¬ 
ences  a  tensile  force.  (Note  a  reversal  of  the 
force  direction.) 

•  The  final  stage  (14  <  r  <  tj)  is  the  vibration 
stage,  during  which  the  plate  exerts  an  oscil¬ 
lating  force  on  the  indentor  where  contact 
may  be  lost  and  regained  many  times. 

Figure  7  shows  the  cumulative  energy  absorption  vs 
time.  Energy  absorbed  increases  rapidly  with  time  in  the 
first  two  stages  and  very  slowly  in  the  third  stage.  During 
the  plate  separation  stage,  a  portion  of  the  energy  returns 
to  the  incident  bar  as  the  indentor  accelerates  into  the  free 
space  in  front  of  it.  This  energy  is  subtracted  from  the 
cumulative  energy,  so  the  energy  curve  dips  at  this  stage. 
Energy  absorbed  in  the  final  stage  is  constant  and  is  the 
energy  required  to  cause  plate  damage. 

Figure  6  also  shows  data  from  two  strain  gauges 
mounted  at  two  positions  on  the  rear  surface.  Strain 
gauge  1  is  mounted  at  the  plate  center  and  gauge  2  near 
the  comer.  The  strain  gauge  data  show  that  plate  bending 
continually  increases  but  remains  small  during  the  first 
four  stages  of  impact.  The  increase  becomes  more  pro¬ 
nounced  in  the  final  vibration  stage  and  results  in  strain 
gauge  failure  near  t  =  ty 

Figure  8  shows  the  laminate’s  energy  absorption 
vs  the  maximum  velocity  of  the  indentor.  The  data 
follow  the  relationship 

£a  =  «  V'  (20) 

where  exponent  m  is  found  to  be  2.55  for  room  tempera¬ 
ture  (21°C)  and  4.53  for  low  temperature  (-30°C).  The 
value  of  the  constant  n  for  21°C  is  3.75  x  10"^  and  for 
-30“C  is  1.01x10"’. 

These  data  are  currently  considered  a  precursor  to  a 
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Figure  7.  Bar  stresses  and  energy  absorption  history. 
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Figure  8.  Influence  of  temperature  and  velocity  on  energy  absorption. 


detailed  analysis  of  temperature’s  influence  on  the  dam¬ 
age  mechanisms  associated  with  composite  plate  im¬ 
pact.  It  is  well  known  that  low  temperature  reduces  the 
fracture  toughness  of  most  materials,  but  no  extensive 
data  are  available  for  composites.  In  research  performed 
by  Nishijima  and  Okada  (1982),  cloth-reinforced  ep¬ 
oxide  resin  composites  were  subjected  to  Charpy  impact 
tests  at  room  and  liquid  nitrogen  temperatures.  These 
tests  demonstrated  that  a  reduction  in  temperature  causes 
a  reduction  in  impact  strength.  Therefore,  for  a  given 


impact  velocity  more  damage  occurs  at  lower  tempera¬ 
tures  and  more  energy  is  absorbed.  The  experimental 
data  in  Figure  8  clearly  show  this  trend  after  a  transition 
velocity  of  about  4.8  m  s“'  (190  in.  s“').  Above  this 
velocity,  more  energy  is  absorbed  by  specimens  held  at 
low  temperature.  Below  this  transition  velocity,  the  low 
temperature  has  an  opposite  effect,  as  discussed  next. 

Greszczuk  (1982)  aixi  Shivakumaret  al.  (1985)  have 
shown  that  in  the  low  velocity  range,  among  non- 
recoverable  energies,  the  Hertzian  contact  force  is  the 


7 


Figure  9.  Effect  of  temperature  on  Hertzian 
contact  energy. 


Figure  JO.  Energy  loss  as  a  function  of 
impact  energy  (after  Sjoblom  et  al.  1988). 


most  dominating  energy  absorption  mechanism.  From 
eq  1 2  this  energy  is  a  function  of  (/:|  +  Substituting 

Vj  =  0.3  emd  E^  =  206.9  GPa  (30  x  10^  psi)  for  a  steel 
indentorintoeq  13, /t,  =  0.9.67  x  10“^.  Referring  to  eq  14. 
^2  contains  the  term  E^  (the  Young’s  modulus  of  the 
composite  plate  in  the  thickness  direction),  which  is 
influenced  at  low  temperature  by  the  change  in  elastic 
modulus  E^  of  the  resin  matrix.  The  empirical  relation 
established  by  Dutta  (1989)  for  matrix  stiffening  with 
lower  temperature  is 

(7-)  =  E„  {To)  -  2.9  X  103  (r  -  To]  (21) 

where  Tq  is  the  room  temperature  and  T  is  the  low 
temperature  in  °C.  Applying  the  above  relationship  to 
the  rule-of-mixtures  equation  (Tsai  and  Hahn  1980) 

E  = - ilMm -  (22) 

VfE^+V^Ef 


E(Em(T) 

£m(7’)  +  Vn,£f 


(23) 


where  and  Vf  are  the  matrix  and  fiber  volume  ratios, 
and  and  Ef  are  the  matrix  and  fiber  modulus  of 
elasticity.  (T)  indicates  the  prevailing  temperature  of  the 
matrix.  Solving  for  E^  and  setting  =  0.31,  ^  is  de¬ 
termined  via  eq  14.  Figure  9  shows  a  plot  of  (/tj  +  ^2)^^ 
vs  temperature.  It  is  clear  from  this  graph  that  as  the 
temperature  is  reduced  the  value  of  E^  (the  Hertzian 
contact  energy ,  which  is  a  direct  function  of  (it ,  +  k2)^) 
decreases.  The  experimental  results  in  Figure  8  confirm 
this. 

In  a  separate  ballistic  experiment  involving  the  pen¬ 
etration  of  identically  configured  30-ply  graphite  epoxy 
laminates  by  12.7  mm  (0.5  in.)  spheres,  Mayer  and 
Altamirano*  observed  energy  absorption  E^  in  the  form 


yields 
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*  Personal  communication  (1990)  from  A.  Mayer,  Air  Force 
Wright  Laboratories,  Dayton,  Ohio. 


80 


Figure  11.  Energy  loss  as  a  function  of  impact 
velocity  as  obtained  from  SJoblom  et  al.  (1988) 
data. 


=  Co  V  ‘-^89  (24) 

where  V  is  the  impact  velocity  and  Cq  is  a  constant. 
Experimental  velocities  ranged  from  500  to  1500  m  s“’ 
(1600-4900  ft  s“*).  Comparing  the  exponents  of  V  from 
the  three  series  of  experiments — 1 .489  from  the  ballistic 
tests,  2.55  from  the  Hopkinson  bar  room-temperature 
tests,  and  4.53  from  the  Hopkinson  bar  low-temperature 
tests — it  is  clear  that  the  damage  mechanisms  and  energy 
partitioning  are  different  at  higher  velocities  than  under 
Hopkinson  bar  low-velocity  impact.  The  lower  value  of 
the  exponent  suggests  a  failure  mechanism  that  progres¬ 
sively  requires  less  energy  as  the  velocity  is  increased. 
This  trend  was  also  validated  by  Sjoblom  et  al.  (1988) 
(Fig.  10).  The  solid  line  in  this  figure  shows  that  energy 
loss  is  a  function  of  V^.  A  reconstruction  of  the  energy 
loss  vs  velocity  (Fig.  1 1 )  clearly  shows  a  decrease  in  the 
energy  absorption  rate  at  higher  velocities. 

CONCLUSIONS 

The  energy  absorption  of  laminated  plates  under 
normal  impact  generated  by  a  Hopkinson  bar  apparatus 
was  exanriined.  Using  wave  propagation  theory,  the 
energy  absorption  agrees  with  theoretical  predictions. 
The  force,  velocity,  and  energy  were  recorded  through¬ 
out  the  impact  event,  which  encompassed  several  hun¬ 
dred  microseconds.  The  experimental  procedure  pro¬ 
vided  a  controlled  condition  enabling  the  full  analysis  of 
events  during  this  period.  Loss  of  contact  between  the 
impactor  and  laminate  occurred  frequently  during  the 
impact  process  due  to  stress  waves  in  the  plate  reflecting 
back  and  forth  through  the  thickness.  The  energy  bal¬ 
ance  equation  provides  a  powerful  tool  for  analyzing  the 
energy  available  to  cause  laminate  damage. 


Low  temperatures  influence  the  energy  absorption 
characteristics  and  control  the  damage  mechanisms. 
Existence  of  a  transition  velocity  was  observed;  above 
this  velocity  energy  absorption  from  an  impact  is  higher 
than  at  room  temperature.  At  low  temperature,  reduced 
impact  strength  causes  a  larger  volume  of  damage. 
Below  the  transition  velocity  (at  low  temperatures)  the 
laminate  absorbs  less  energy,  which  can  be  approxi¬ 
mated  by  Hertzian  contact  theory. 

When  comparing  data  generated  via  the  Hopkinson 
bar  apparatus  to  data  obtained  via  high  velocity  ballistic 
impact,  the  failure  mechanism  is  seen  to  be  a  function  of 
impact  velocity.  As  the  impact  velocity  increases,  the 
fracture  requires  less  energy. 
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APPENDIX  A:  WAVE  PROPAGATION  IN  A  HOPKINSON  BAR  APPARATUS 

The  incident,  reflected,  and  transmitted  energy  are  denoted  by  and  E-j-  respectively.  The 
governing  equations  for  stress  a{t),  velocity  V(r),  and  displacement  u(t)  are 

(T(f)  =  (p)(c/g)V(/)  (Al) 

V  =  duldt 

i 

tt(f)  =  (c/£)J  a{t^dt 
t=o 

where  E  =  Young’s  modulus, 
t  =  time, 
p  =  density, 

V{t)  =  velocity, 

c  =  the  wave  speed  in  the  material,  and 
g  =  the  acceleration  due  to  gravity. 

Further,  the  energy  of  the  wave  in  a  one-dimensional  bar  is 

t 

£  =  4c  a\t)dt  (A4) 

E] 

t=0 

where  A  is  the  cross-sectional  area  of  the  circular  Hopkinson  bar.  Displacement  of  the  indentor  up) 
can  be  obtained  from  the  incident,  reflected,  and  transmitted  stresses  (0|(r),  and  OjfO. 
respectively): 

I 

Ml  (r)  =  (c/E)  j  [oi  (f)  -  Or  (f)  -  07(1)]  dt  (A5) 

/=o 

The  force  on  the  indentor  is 

E(r)  =  (A)[oi(r) +aR(/)]  (A6) 

Energy  produced  by  the  indentor  (E  )  can  be  computed  from  an  integration  of  the  product  of  force  and 
displacement  over  the  duration  of  the  wave: 

I 

Ep  =  (Ac/E)J  j[ai(t)  +  aR(t)]-[0,(t)-OR{t)- OT(t)]j  dr  (A7) 

f=0 

The  energy  of  the  incident,  reflected,  and  transmitted  waves  can  then  be  computed  from 

I 

(E, ,  Er  ,  Et)  =  {Ac IE)  I  [o,(t)2  ,  OR(t)2,  OT(t)2]  dt  (A8) 

/=0 


(A2) 

(A3) 
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The  total  energy  absorbed  w  ithin  the  system  (i.e.  that  absorbed  during  damage  generation  and  minor 
losses  within  the  fixture)  is 

E^  =  El  -  {Eji,Er),  (A9) 

which  can  be  rewritten  as 

t 

=  {Ac IE)  j  ja,  (f)2  -  [or  (f)2  +  OT(r)2]  j  dt  (AlO) 

(=0 


The  Hopkinson  bars  used  in  this  test  series  were  steel  with  the  following  parameters; 

A  =  1 140  mm^  (1.7675  in^) 
c  =  5017  ms-'  (197531  in.s-') 

E  =  193.05  GPa  (28  X 10^  psi) 

so  that  Ac/E  =  0.02963  m^  GPa-'  s''  (0.012467  in.^  lb”'  s”'). 
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